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Jerome MartinQ Carlo SchimdQ and Jean-Philippe UzaiH 

Institut d'Astrophysique de Paris, QReCO-CNRS, UMR 7095, 
Universite Pierre et Marie Curie, 98bis boulevard Arago, 75014 Paris, France 
(Dated: February 5, 2008) 

In scalar-tensor theories of gravity, the equation of state of dark energy, w, can become smaller 
than —1 without violating any energy condition. The value of w today is tied to the level of 
deviations from general relativity which, in turn, is constrained by solar system and pulsars timing 
experiments. The conditions on these local constraints for w to be significantly less than —1 are 
established. It is demonstrated that this requires to consider theories that differ from the Jordan- 
Fierz-Brans-Dicke theory and that involve either a steep coupling function or a steep potential. It 
is also shown how a robust measurement of w could probe scalar-tensor theories. 
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Various observations indicate that the expansion of our 
Universe is presently accelerated pj . While still debated, 
this conclusion appears to be more and more robust and, 
as a consequence, the discussion has now mainly shifted 
to explaining the cause of this acceleration. The property 
of the effective equation of state inferred from the obser- 
vations, w, is a key issue in this investigation. In particu- 
lar, showing that w ^ — 1 and/or dw/dz ^ (z being the 
redshift) would exclude a cosmological constant, proba- 
bly the most natural candidate and, hence, would have 
drastic implications for fundamental physics. Recently, 
various observations have pointed towards the conclu- 
sion that w < — 1 |2| • Although far from being settled, 
this would also have important consequences since such 
an equation of state cannot be achieved by quintessence 
models for which — 1 < w < 1. "Phantom models" Q, 
consisting in a scalar field with a minus sign in front of the 
kinetic term, are very often advocated in order to explain 
W < — 1. These models are plagued by various theoret- 
ical problems such as, for instance, their stability when 
interactions with other fields are taken into account. 

However, another route can be investigated since theo- 
ries where the gravity sector is modified, i.e. where grav- 
ity is no longer described by general relativity (GR), can 
also entertain w < — 1 even if the matter sector is de- 
scribed in a standard fashion. The prototype of such a 
theory is a scalar-tensor theory of gravity which is both 
well-defined and well-motivated H as they arise as the 
low energy limit of string theory. It is worth recalling 
that they are equivalent to theories where the gravita- 
tional action is given by an arbitrary function of the Ricci 
scalar and also encompass the Jordan-Fierz-Brans-Dicke 
(JFBD) theory as a particular case. As already men- 
tioned, having 1 + w negative is linked to the fact that 
the gravity sector is modified but such modifications are 
strongly constrained by solar system and pulsars timing 
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experiments. Therefore, one can hope to use these local 
tests to track the real nature of the dark sector. In this 
letter, we investigate these issues in general scalar-tensor 
theories. 

In the Jordan frame, i.e. in the (physical) frame where 
the experimental data have their usual interpretation, 
scalar-tensor theories are described by the action |(j 

-2U{<p)]+S m [9 m , giu/ ], (1) 

which depends on three arbitrary functions F (tp), Z (tp) 
and the potential U (tp) , only two of which are indepen- 
dent. G* is a constant, different from the gravitational 
constant measured in a Cavendish experiment, G cav , as 
will be discussed in more details below. In the following, 
the matter action, S m \$f m , guv] i describes a pressure- less 
perfect fluid. 

In the particular case of a Friedmann-Lemaitre- 
Roberston- Walker (FLRW) Universe, choosing Z (tp) = 
1, the field equations reduce to 

3 (h 2 + ;|) = 8 „ G ,^ + £- 3 4 + £, (2) 

fa 9 K.\ tp 2 F F U , . 

-[2-+H 2 + ^r) = — + - + 2H ,3 

\ a a 2 J 2F F F F' v; 

with H = a I a, a being the scale factor and K, the 
curvature of the spatial sections. The energy density 
of matter scales as a -3 . These equations should be 
compared to 3(H 2 + JC/a 2 ) — 87rG cav (p m + p DE ) and 
— (2a/ a + H 2 + K./a 2 ) = 8nG cav p DE , since the equa- 
tion of state of dark energy is experimentally inferred 
from the expansion history of the Universe by using the 
standard Friedmann equations of GR. It follows that 
3f2 DE u> = — 1 + ^ic + 2g where q = —aa/(a) 2 is the 
acceleration parameter, tt/c = —K./a 2 H 2 and that the 
dark energy density parameter is defined as Sl DE (z) = 
H 2 /Hl - ng^l + zf - Q° K (1 + z) 2 , where here and in the 
rest of this letter the subscript "0" denotes the present 
day value of the corresponding quantity. Clearly, iden- 
tifying p DE and p DE from Eqs. J2J and © leads to the 
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conclusion that, in scalar-tensor theories, w needs not to 
be positive definite, depending on the choice of F. This 
was already stressed by various works Q and worked out 
in the particular case of JFBD theories in Refs. 

The present values of F and its derivatives are con- 
strained by the solar system and pulsars experiments. 
In order to reveal the link between w and these con- 
straints, it is more convenient to work in the Einstein 
frame. The Einstein frame metric, g* v (all Einstein frame 
quantities will be denoted with a star), is related to the 
physical metric through the conformal transformation 
9% = F (f)9^- Setting F(<p) = A -2 (</>*), the scale 
factors and cosmic times in both frames are related by 
a = Aa* and dt — Adt*, so that the Hubble parameters 
are linked by AH — H* (1 + ap'^) with H* = dlna*/di* 
and where a prime denotes a derivative with respect to 
the Einstein frame number of e-folds, N* = ln(a*/a*o)- 
Deviations from GR are usually described by two pa- 
rameters, a and /3, which are the first and second field 
derivative of the coupling function A(ip*), namely 
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In this framework, GR is characterized by a = (3 = 
while the JFBD theory corresponds to a constant a 2 = 
(2u BB + 3) _1 and (3 — 0. As mentioned before, in 
scalar-tensor theories, the Newton constant obtained in 



a Cavendish experiment differs from G* and is given by 

G cav = G*A 2 o(l + al) . (5) 

The parameters ao and (3q are constrained by various ex- 
periments. If we define 7 = 7 PPN — 1 and (3 = f3 rpN — 1, 
where 7 PPN and /3 PPN are the post-Newtonian parame- 
ters 0, related to ao and /3 by 7 = —2oq/{1 + a 2 ,) and 
2(3 = (3 al/(l + al) 2 , then the perihelion shift of Mer- 
cury implies 2^ — /3\ < 3 x 10~ 3 while the Lunar Laser 
Ranging experiment sets 4/3 — 7 = —(0.7 ± 1) x 10~ 3 . 
On the other hand, a bound on 7 alone is set from the 
time delay variation to the Cassini spacecraft near solar 
conjunction, namely 7 = (2.1 ± 2.3) x 10~ 5 , see Ref. 01 
for a review. We conclude that 
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where the lower bound on (3$ arises from pulsar timing 
experiments |l0j . Note, however, that we cannot con- 
sider arbitrarily large values of (3q since then the post- 
Newtonian approximation scheme would breakdown. In 
this case, the upper constraints should be reanalyzed. We 
thus loosely assume that (3q < 100. 

Working out the Friedmann equations in Einstein 
frame 0, we obtain that the equation of state is given 

by 
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where n m = 8TrG caiV p m /3H 2 so that Q DE +fl m + Q,/c = 1. 
In the limit of a minimally coupled scalar field, it reduces 
to the standard relation 1 + w = 2p' 2 /S£l DE . 

At this point, it is of utmost importance to stress that 
the value of p'^ is not free. Indeed, in the Einstein frame, 
the Friedmann equation reads Q 
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+ 87rG*p m » + 2V(tp*) . 



(8) 



where V(pA — U/(2F 2 ), and the positivity of the energy 
density of matter implies that \<p'A < \/3, as long as 
flic ri m (in the following, we assume that flic — 0). 

The expression Q for w is completely general but very 
intricate and hence not so illuminating. However, taking 
into account that ao has to be small and that p'^ has to 
be bounded by v3 and, therefore, that ao^io has to be 
small as well, the present day value of the equation of 
state simplifies considerably and reduces to 

3f!° E (l + w ) ~ 2(1 - /3oV?o ~ W*V (9) 

This formula turns out to be our main result. The contri- 
bution of Pq arises from the term F/F in the right hand 



side of Eq. ©. In fact, Eq. JjJ) shows that wq < — 1 
is always possible provided (3q > 1 (in this case, even if 
the slow-roll approximation is satisfied, i.e. <p" <C p'*o) 
and/or aip'l positive and large compared to ip' 2 . Both 
regimes cannot be reached in the case of JFBD theories 
(except if the time variation of G cav is large, see below) 
and exist even in the limit ao — » so that all local tests 
can be satisfied. 

The amplitude of Wq depends on the value of ip' r0 and 
V>"q. Independently of any dynamics, these two quantities 
are constrained by the bounds on the time variation of 
the gravitational constant ^l| . Defining 
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(10) 



the parameter ao is bounded by |o"o| < 5.86 x 10~ 2 /i -1 . 
There is no stringent bounds on £0 but we can estimate 
that, since ao has been "measured" during a period of 
about 20 years, we have |£oi?ol ^ |G/G| /(20 yr). This 
implies that |£ | < 5 x lO 8 /!" 1 ^ - 2.5 x 10 7 h~ 2 . Using 
Eq. JSJ, one can then express p>'^ and p'l in terms of the 
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parameters a and £. As long as ^ — (1 + a 2 ), a case while the second derivative reads 
we shall discuss later, one arrives at 

.of a\ 1 
2a \ 1 + cr 2 / 
I 



2a[l + /3/(l + a 2 )] 



-(i 



1 + a 2 , 



a + 



i 



l + a 2 +/3 \dtp 



d/3 2a/3 



1 



1 



(12) 



Interestingly enough, we see that <p" , and hence Wq, de- 
pends on £o and on the derivative of the function. 
The latter is not constrained and we will assume that 
d/3/d<y9* < 0(100) for the same reason as for 0. The 
above formula is rather complicated but, given the pre- 
vious constraints, we simply have 
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In particular, the unknown term dft/dip* does not appear 
in this approximation. 

As mentioned above, the previous considerations are 
independent from the dynamics. However, from a model 
building point of view, tp" and (/?'„ are not independent 
once the potential V((f*) is chosen. They are related 
through the Klein- Gordon equation which reads 

y + i } <p* + ( 2 + x h = -( ax + , (w) 

with X = n m /il u: where p v = 2U/(16nG*) and a v 
( 1 In V/dip* . When the kinetic energy of the scalar field is 
negligible, we have Q v ~ f2 DE . 

We now come back to our master equation © and an- 
alyze the two regimes where 1 + wq can become negative 
and even large (in absolute value). The first regime is 
the natural extension of quintessence models to scalar- 
tensor theories |l2j and corresponds to the situation in 
which the field is decelerating, i.e. ip" <C (p'* , so that 
O!o<p" is negligible in Eq. (JSJ. Clearly, this requires 
fto > 1. Such a regime, that cannot be reached in a 
JFBD model, has the advantage to exhibit an attraction 
mechanism toward GR |13| so that ao can dynamically 
be made small. Since, in this situation, 0q is not close 
to — 1, Eq. (|TT)l implies that ip'^ ~ a /[2a (l + O )] and 
Eq. (JT3J) leads to £ ~ /Vo/pa^l + ft )}. On the other 
hand, from the Klein-Gordon equation, one obtains that 
</?*o ~ -(ao^o + avo)/(2 + ^o) an d, in order to fulfill the 
condition ip'^Q < y/3, one must have |a v o| ^ 4.2. The two 



expressions for p' r0 can be used to infer what fto is and, 
if we insert the result in Eq. ©, one gets 
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(15) 



in the most interesting limit where a v o 3> aoXo since, 
otherwise, |1 + 10o| ~ 0(oq). They are two ways to inter- 
pret Eq. (IS) . Either it gives wq in terms of ao and a v o> 
assuming that ao is known. Or it provides, for fixed ao 
and a v o, the minimum value that wo can reach assum- 
ing, as is the case now, that only an upper bound on ao is 
available. Fig. ^ depicts the above-mentioned minimum 
value as a function of (ao, avo). 

At this point several remarks are in order. Firstly, 
it is unclear whether the above situation can be real- 
ized without some fine-tuning in a realistic model, where 
both the coupling function and the potential are running 
away in order to have attraction towards GR and insen- 
sitivity to the initial conditions [l4|]. Let us illustrate 
this point on the following particular example. Consider 
the case where a cx e~ x,p * so that fto = — aoA. Since 
we need ao « 1 and, at the same time, fto > 1, this 
means that A ^> 1 which may be considered as unnatu- 
ral. Secondly, one can reverse the logics and study what 
a detection of wo < —1 would imply on scalar-tensor 
theories (assuming, of course, that a slow-rolling 92* is 
causing the acceleration of the expansion). Besides the 
fact that the condition 0o > 1 would drastically improve 
the current limit on 0o and, in particular, exclude GR, 
it is also interesting to remark that this would link wq 
to the time variation of the Newton constant through 
the expressions Oq — 6 ^d E ( 1 + w o)(l + /3o) 2 «o/(l - fto) 
and £0 ~ 3n° DE o (fto + 1)(wq + 1)/(1 - fto) < 6 x 10 2 , a 
limit sharper than the bound set by local experiments. 
Finally, it is interesting to notice that wq < ~ 1 can be 
obtained even if the potential energy is negligible (in GR 
this would correspond to wo ~ 1). In this case, one 

can show that <^ ~ (S^J 1 / 2 and ip% 3Q°y* /2 

so that aQ(p" is indeed negligible in Eq. (J5J). Then, 
the equation of state can easily be obtained and reads 
1 + wo ~ 2(1 — fto). Again fto is bounded by the con- 
straint on ao, l + 0o < coA^o^/lSS!^), so that a large 
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FIG. 1: Minimum value of wo as a function of (qo,q v o) for 
<7 < 1(T 3 . 



value of wq requires <7n <C an • 

The second regime corresponds to the case where oup'l 
dominates in Eq. (J5J). As can be seen in Eq. (|13fl . this 
is possible if £n is large and positive and/or if /3n ~ — 1. 
Strictly speaking, one should in fact consider the limit 
[3 — > — (1 + an) corresponding, when V — 0, to the 
Barker theory 15] in which A = cos <p* and G cav con- 
stant (er = £ = 0, whatever the value of ip'^ and ip"). 
Since a = — tan tp* , the cosmological evolution drives the 
theory away from GR unless a potential keeps ip* close 
to until the last e-folds inducing a large variation of 
aip" ~ Oil) in order to have 1 + w < 0. Such a model 
seems very contrived unless the potential exhibits a slope 
discontinuity recently. 

Let us now assume that £n ~ since, from dimen- 
sional analysis, one expects ^ ~ ff o ~ ^ x 10~ 3 . In 
the limit f3 n — > — 1 — eft, Eq. (fHl)l leads to ap'l ~ 
2ao/3o¥'*o/(l + 00 + «o) an d wo diverges when /3n — > 
— 1 — Q!q. However, the constraint <p' % < V3 implies that 



1 + /?o + ccq <; <To/(2\/3ao) and, hence, the smallest value 
that can be obtained is (1 + iun) m i„ ~ — 8v3ao/ (^ e ^o)- 
In particular, wq <C — 1 is perfectly possible in this regime 
if Cn -C Uq. Finally, let us see what this regime implies 
in terms of model building. The Klein- Gordon equation 
implies that 2<p" /3 ~ — a v n/(2 + Xn) so that the slope of 
the potential must be very large, a v o 3> • Moreover 
since, in this regime, /3n < the theory is driven away 
from GR and it follows that the potential must be tuned 
in order to prevent this drift. 

As a conclusion, let us summarize our main findings. 
Firstly, we have confirmed that 1 + u>n is not positive 
definite in scalar-tensor theories even if all the matter 
energy conditions are satisfied. Secondly, we have estab- 
lished under which conditions 1+wq can become negative 
given the local constraints coming from solar system and 
pulsars measurements (in the case of chameleon mod- 
els [iff , it has been argued that an can be of order unity 
today, a case where large negative values of 1 + wq can 
be achieved more easily). We have shown that getting 
a non negligible deviation from —1 necessarily implies a 
non vanishing /3n (except if £n is large), a situation that 
cannot be reached in the JFBD case. Thirdly, in terms 
of model building, we have demonstrated that this cor- 
responds either to /3q > 1, a situation where the scalar 
field is slow-rolling today and the coupling constant is 
very steep or to (3q ~ — 1 (or £ large) where the slope of 
the potential is very large. Fourthly, we have also shown 
how a measurement of wo < —1 could improve the local 
constraints on the deviations from GR. This highlights 
the complementarity [TtI ] between cosmological and lo- 
cal tests of gravity. Finally, let us remark that we have 
only taken into account the local constraints. A next 
step would be to reconstruct the redshift evolution of the 
models and to show that they are not pathological [||. 
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